AEI-058 



|hep-th/9803TT7 



Monte Carlo Approach to M-Theory 



Werner Krauth * 
CNRS-Laboratoire de Physique Statistique, Ecole Normale Superieure 

24, rue Lhomond 
F- 75231 Paris Cedex 05, France 



Hermann Nicolai ' § and Matthias Staudacher * § 
Albert-Einstein-Institut, Max- Planck- Institut fur Gravitationsphysik 

Schlaatzweg 1 
D- 14473 Potsdam, Germany 



Abstract 

We discuss supersymmetric Yang-Mills theory dimensionally reduced to zero 
dimensions and evaluate the SU(2) and SU(3) partition functions by Monte 
Carlo methods. The exactly known SU(2) results are reproduced to very 
high precision. Our calculations for SU(3) agree closely with an extension 
of a conjecture due to Green and Gutperle concerning the exact value of the 
SU(N) partition functions. 



*krauth@physique. ens.fr 

t nicolai@aei-potsdam.mpg.de 

* matt hias@ aei-p otsdam . mpg . de 

Supported in part by EU Contract FMRX-CT96-0012 



1 



I. THE MODEL 



Recently there has been renewed interest in dimensionally reduced supersymmetric Yang- 
Mills theories. It has been demonstrated that reductions of the ten-dimensional theory to 
p+1 dimensions are relevant to the description of p-dimensional extended objects (D-branes) 
in string theory, and therefore, potentially, to quantum gravity. Furthermore, various such 
reductions are currently being investigated in attempts to find a theory, termed M-theory, 
which, it is hoped, non-perturbatively encompasses all known superstring models as well 
as D = 11 supergravity. The most extreme such reduction, to zero dimensions, is believed 
to describe string instantons, i.e. configurations of open superstrings whose ends are fixed 
in space-time. The functional integral of Yang-Mills theory becomes an ordinary multi- 
dimensional integral, which turns out to be well-defined and finite if supersymmetry is 
present. The same integral appears in attempts to rigorously define and calculate the Witten 
index of quantum- mechanical supersymmetric gauge theory Lastly, it is the cornerstone 
of the so-called IKKT model |3|] which is claimed to yield a non-perturbative definition of 
superstring theory. Explicitly, the integral reads, for gauge group SU(N), 




(1) 

We do not include a coupling constant into this definition, since it could be immediately 
scaled out in a trivial fashion. The matrices in the exponent in eq. ([!]) are in the fundamental 
representation of SU(N), i.e. = X a Ta, ^ a = *& a Ta, where the SU(N) generators T4 
are hermitean and normalized such that Tr T A T B = \5 AB '. The symmetric M x jV matrices 
are related to the standard SO(l,D — 1) gamma matrices by V 1 = Cj^, where C is 
the charge conjugation matrix. The model is supersymmetric in dimensions D = 3,4,6, 10, 
where the degree J\f of (real) supersymmetry is, respectively, Af = 2(D — 2) = 2,4,8,16, 
with the supersymmetry variations 

6X li = H'f* 5V = - l -[X„X u ]Yre. (2) 

The most interesting case, in view of the above applications, is the maximally supersym- 
metric integral with D = 10. 

We note that all calculations are performed with a Euclidean signature. The Wick 
rotation to Euclidean space is simply accomplished by the substitution Xq — * iXo- This 
leads to a positive definite bosonic "action" in the exponent, which is a necessary (but not 
sufficient) prerequisite for the integral to be well-defined. When switching from Minkowskian 
to a Euclidean signature, one in principle faces the problem of fermion doubling because 
e.g. the Majorana-Weyl spinors required for maximally supersymmetric Yang-Mills theory 
exist only for signature (1,9). However, this problem is spurious and completely resolved 
by working only with the Euclidean Majorana field and never making use of its complex 
conjugate field [|]. This is consistent because the relevant requirement in Euclidean quantum 
field theory is not hermiticity, but rather Osterwalder-Schrader reflection positivity. 

Integrating out the real Grassmann variables ty A , one obtains 
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r N 2 -l D fiX A pi -, 

z d,n= / II Il^= ex P ^[X^X^X^Xv] V D , N (X), (3) 
J a=i M =i v2vr U J 

where Pd.tv is a homogeneous polynomial of degree k = (D — 2)(N 2 — 1) = ^Af(N 2 — 1) in 
the variables X A . Its crucial property is invariance under SO(D) x SU(N). This polynomial 
is just the Pfaffian of the 2k x 2k antisymmetric matrix 

{M D , N )t* = ^f ABC K,XC, (4) 

which can be directly read off from eq.(^). Here f ABC = — 2iTr T A [T B , T c ] are the SU(N) 
structure constants. The Pfaffian is explicitly given by 

V D , N = PiM D>N = Ei-^M^ • • • M2y U2k - (5) 

Here a runs over all permutations of the 2k double indices a = (aA), and (— 1)°" is the sign 
of the permutation. Although the matrices M. o,n are complex in general, the Pfaffians are 
real as a consequence of the fact that we started from a hermitean action (in Minkowski sig- 
nature). However, they are not necessarily positive. The explicit form of A4d,n depends on 
the particular representation of the Clifford algebra and the corresponding gamma matrices 
T^. For defmiteness, let us specify a convenient choice for the four dimensions of interest. 
Denote by the adjoint representation of SU(N), i.e. X AB = f ABC X^. 

•D = 3: Here the Pfaffian is determined from the 2(N 2 — 1) x 2(N 2 — 1) matrix 

x«(a-)^( x ^ x ; x ; X2 ) (6) 

•D = 4: In this case it is possible to reduce to a complex two-dimensional representation 
instead of the real four- dimensional representation. This allows us to rewrite the Pfaffian 
as an ordinary determinant: there exists a 2(N 2 — 1) x 2{N 2 — 1) matrix -M 4A r such that 
Pf_M 4i 7v = detM' 4:N , with 

•D = 6: Again it is possible to replace the real eight-dimensional by a complex four- 
dimensional representation: we have PfA^v = detA4' 6N with the 4(iV 2 — 1) x 4(iV 2 — 1) 
matrix A4' 6 N 



M' N (X) 



X 6 + iX 5 iX 4 + 01X1 + <x 2 X 2 + (T3X3 

iX 4 - cr 1 X 1 - <x 2 X 2 - 03X3 X 6 - iX 5 



where the <jj are the standard Pauli matrices. 

•D = 10: In this case we use the decomposition 



M 10 , N (X) - ( zXr _ px Xiq _ zXg ] (9) 
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for the 16(iV 2 — 1) x 16(iV 2 — 1) matrix M.io,n- The indices i, j, k run from 1 to 7, and the 
seven 8x8 matrices P satisfy the 5*0(7) Clifford algebra {r j ,r jf } = 25*7 A particularly 
nice representation is in terms of the octonionic structure constants (again with i, j, k = 
1,..,7) 

(r*) ifc = ici ifc (P) 88 = o nm 
(T%= -(T% = i<5j liUj 

These structure constants are completely antisymmetric; the non-vanishing ones are fully 
specified by: 

Cl24 = C235 = C346 = C457 = C56I = C672 = C713 = 1 (11) 

The P are obviously hermitean and purely imaginary. Note that, unlike for D = 4 and 
D — 6, it is generally not possible here to reduce the Pfaffian to a determinant. The only 
exceptions are N = 2 and N = 3, since in this case the number of dimensions (=10) exceeds 
the dimension of the gauge group, and we can use 50(10) invariance to set Xg = X w = 0. 
Then the Pfaffian becomes 

Vio,n(X)\ X9=Xio=0 = det(*X 8 + PX.) (12) 

This allows to integrate, for SU (3) and D = 10, over only eight instead of ten matrices, an 
observation that will prove useful in the numerical calculations below. After proper inclusion 
of the Fadeev-Popov determinant associated with this reduction, the integral becomes 



/8 n]X A Tl 8 1 n2 

IJ ^=exp -Tr [X a ,X b ][X a} X h ] det^c^f detM'(X) (13) 



A=1 V2tt 



a,b=l 



where M'{X) = *X 8 + PXi. 

For general D and N, the Pfaffian polynomials Vd,n have a very complicated structure 
about which little is known. These complications are not yet fully apparent for SU (2), where 
the Pfaffians are still comparatively simple |],[2|,[5[]: 

Vs, 2 (X) = \e ABC e, vp X^X c p = (14) 

= ~ie lwp TrX ll [X v ,X p ] (15) 

p 4>2 (X) = \^^XiXi'X^XfXf = (16) 

= ^Tr [X„ X v \ [X u , X p ] [X p , X,} (17) 

Vw(X) = (^Ti[X p ,X u ][X u ,X p ][X p ,X,]) 2 (18) 

V W> 2{X) = (^Tr [X p ,X u ][X u , X p ][X p , X p ])' (19) 

Remember that the Greek indices in the preceding four equations are summed from 1 to, 
respectively, 3,4,6, 10. Already for 577(3) and D = 3, however, more ingenuity is required; 
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using symbolic computation, we have found the result (the d ABC are the SU(3) symmetric 
structure constants) 



P 3)3 (X) = (3.fABE d CDE f A'B>E> d C>D'Ei_ 

3 fABC jDA'E fEB'E' jE'C'D'\ v A v A' v B v B' v C v C' y D v D' 
— 10/ Cl J d I A p A p Ji v Ji y A p A p A CT A CT 



or equivalently, 



I TrfX^X.KX^Xj TrX^X^X^Xj + 
6 ^ v [X v 



■ w, -X/j] TrX^fjXj,, X CT }, {X p , X CT }] 



(20) 



(21) 



For higher dimensions and yet larger values of X, we expect the Pfaffians to increase 
rapidly in complexity. Let us nevertheless collect the following general relations whose 
validity follows upon inspection of the explicit form of the matrices M.d,n given above. 
Namely, we have, for any X, 



^10,Af(^)lx 7 =...=X 10 =0 
^6,iVpOlx5=X 6 =0 

"P-MvPO lx 4 =o 



(7\at(X)) 2 
(7VP0) 2 

(?VP0) 2 



(22) 



Of course, these relations are consistent with the explicit SU{2) results quoted above. We 
emphasize, however, that relations of this type do not hold for generic configurations, as we 
have verified numerically. In addition, for X > 3 and D > 6, we have found configurations 
X M whose associated Paffians are negative. We also observe that Vd,n{X) vanishes whenever 
the matrices X M are restricted to lie in the regularly embedded su(N — 1) subalgebra of 
su(N). Finally, Vd,n{X) vanishes if only two X^'s are different from zero. 

Returning to the integral (3), the first question to ask is whether or not it exists. The 
following naive argument seems to indicate that Zd,n is divergent: Diagonalizing one of the 
matrices X M and restricting to the subspace where the off-diagonal elements of the other 
D — 1 matrices are zero as well, we see that the diagonal components of all D matrices do 
not appear in the action anymore. One might therefore suspect that the fluctuations along 
the diagonals cause the integral to diverge. Explicit calculations disprove this reasoning: it 
is shown in that for the gauge group SU (2) the exact values of the integral (j|) areQ 



Z 



D,2 



57T X < 



'0 


D 


= 3 


1 

4 

< 


D 


= 4 


1 
4 


D 


= 6 


5 

, 4 


D 


= 10 



(23) 



1 While all the technical ingredients for this calculation can already be found in |fj], the final 
result for the index obtained there differs from the correct result in eq.(]23|) by a factor of 4. Also, 
that calculation did not take into account the deficit term of || needed for the interpretation as a 
Witten index. 
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It is an instructive exercise to repeat the SU(2) calculations for the case of no supersymmetry 
(A/" = 0). Omitting the Pfaffian, one obtains 



r 3 D dX A rl 
/ II n^=exp -Tr[X^X v ][X^X v ] 

J A_i VZ7T l£ 



A=l ju=l 



oo D < 4 

(24) 



_3 D _ i r (£)r(V)r(V) 



r(f)r(^i)r(^2) 



One concludes that despite the valleys of the potential, entropic effects of the measure 
can become strong enough to overwhelm the possible divergences. Therefore we conjecture 
the integral Z^ ^ to be finite for all N, even though no rigorous proof is available for 
N > 3 so far. There are two sources of such entropic effects. One is the increase in 
dimension with larger D and N: In the above example, the integral becomes convergent 
if only a sufficient (here: five) number of random matrices is added. Furthermore, we 
generally expect the integrals to become better behaved with larger N since the relative 
weight of the valley configurations with regard to the generic configurations decreases like 
N^ 1 . Secondly, convergence is improved by supersymmetry because the Pfaffian damps the 
contribution of the valleys: It is easy to see that the Pfaffian vanishes if all D matrices 
commute simultaneously. 

Coming back to the supersymmetric case, we would like to find the generalization of 
eq. (|23"D to higher rank gauge groups. Unfortunately the methods of are very special 
to SU(2) and cannot be generalized to N > 2 in any obvious way. On the other hand, 
exploiting the suspected relationship of eq.(||) to superstring instantons, Green and Gutperle 
made a bold conjecture as to the value of the D = 10 integral. Noticing, for N = 2, 
a numerical coincidence between the result ( p3[) and a term in their calculation of the D- 
instanton effective action, they suggested 

Zd,n ~ E \ for D = 10 ( 25 ) 

Here the sum runs over all divisors of N. If N is prime, the sum is simply 1 + j^. Now 
it is clear that such a conjecture is meaningless unless one very carefully specifies all the 
normalizations in eq.(|3|)! Indeed, both the action and the polynomial Vd,n are homogeneous 
functions of the X£, and therefore @ is a pure number up to rescalings: It is impossible 
to introduce any non-trivial coupling constants. Fortunately, a natural normalization is 
suggested by the index calculations of 0. An iV-dependent "group factor" Tn has to be 
included, which gives precisely the extra factor JF 2 = \/8ti in (p3|). We can then extend the 
conjecture of to also include the remaining dimensions D = 3,4, 6, and propose^ 



2 This proposal slightly differs, for D = 4, 6 and non-prime values of N, from an earlier version 
of the present article. We thank I.K. Kostov for alerting us to the fact that the present form is 
indeed more natural. 
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T N x < 



n 
\j 


D 


— 3 


1 

N 2 


D 


= 4 


1 

N 2 


D 


= 6 


V , J- 

£—/m\N m 2 


D 


= 10 



(26) 



We note that the result Z^^ = for even N is a trivial consequence of the reflection 
antisymmetry — > — X M ; on the other hand, as our formulas (|2"U|) , (|2"T|) for the D = 3, N = 3 
Pfaffian shows, the result (if true!) is all but trivial for odd N. 

It seems natural to assume that Tn is independent of the dimension. Extending the 
assumptions going into the index calculation of 0] to arbitrary iV > 2 we find, with our 
normalizations of the integral @, 

N(N+1) Y-l 
~- 2 2 7T 2 

= - 7=-^ ., (27) 



It should be stressed that our extended conjecture (pB]), if true, furnishes an esaci analytic 
solution of the IKKT matrix model ||. Indeed, the partition function of the IKKT model 
is defined by summing up the D = 10 SU (N) partition funtions, and we thus propose 

oo 

ZjkMP) = E Z io,n e~ pN (28) 

with i?io,7v given as in eqs. (p6|) , (p7|) , and i?io,o = ^10,1 = 1- The weakest part in our proposal 
is the form ( P?[ ) of (see also our results below): Firstly, we might still have missed an 
iV-dependent factor in eq.(|27|), and secondly, the correct definition of the sum over iV might 
involve a yet unknown iV-dependent constant 

OO 1 

2ikkt(/3) = E C N E — 2 e~ m (29) 

N=Q m \N " l 



It is interesting to observe that for D = 10 the proposed form (|2q ) of Zo,n is non- analytic 
in iV: For any real number c with 1 < c < ^- we can find a sequence {N w } C IN such that 
lim^oo Z)m|Ar„ ^2 = c. Such a result, if true, would give an explicit example of a case where 
even the assumption of the existence of a large N limit is simply wrong. On the other hand, 
it is intriguing for the future that summing over all values of N, as proposed for the IKKT 
model, clearly reinstates analyticity (as long as Cn is analytic), since we can rewrite eq.(p9|) 
quite generally as 

OO 1 oo 

Z 1KKT ((3) = — E C mk e- pmh (30) 

Clearly these conjectures are very much in need of direct evidence from the actual integral 
(§). Analytical methods to treat this type of integral for N > 2 are not available for the 
moment, and we have therefore concentrated on a numerical evaluation of the SU(2) and 
ST/ (3) partition functions using Monte Carlo methods, which, as we will show, need to be 
quite intricate in order to cope with the large variations of the integrand, provoked by the 
presence of the valleys. 
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II. MONTE CARLO APPROACH 



In the following, to simplify notations, we write the integral eq.(H) as 



Z 



D,N 



/ dx z D)N (x) (31) 

J \F 



where d = D(N 2 — 1) is the total dimension of the integral, with x = (x\,X2, ■ ■ ■ , Xd)- The 
volume V, for the moment, is V = IR d . For a Monte Carlo evaluation of widely varying 
integrals (cf |7J]), we generally attempt to split off a positive importance sampling function 
ir(x): zd,n(x) = n(x) x zd,n{ x ) such that 

J dx Zd,n{x) — J dxir(x) z D:N (x). (32) 

A judicious choice of it will permit us to preferentially generate points in the valleys. This 
will much reduce the statistical fluctuations: Suppose we are able to draw M independent 
random configurations Xi, i — 1, . . . , M distributed according to the probability p(x) ~ tt(x) 
with Vt, = j v dxn(x) . We then have 



1 \S 2 Z DN 
2d,n = K x — zdA*) + °( ITT ) (33) 



with 



M ^ M 



5 2 Z D , N = V n 2 x- ^d,n(^)) 2 - V* x E ZdA*))' (34) 

i=l i=l 

The unavoidable statistical error will therefore depend on the second moment 

dxTr(x) z 2 D N (x) (35) 



v 

In principle, the Monte Carlo computation is possible only with choices of ir for which the 
second moment (^) exists. It should also be realized that, as is evident from eq.(^), we 
need to know V^, the integral of 7T. If the latter is unknown, we are usually only able to 
evaluate expressions of the type 

Zp tN _ Jdxw(x) z DiN (x) 
V n J dxir(x) 

We summarize these last remarks by stating that every approximation n(x) to the integrand 
zd n in eq. (§) can serve as a useful starting point for very precise Monte Carlo evaluations 
of the integral - for arbitrarily high dimension d - provided that J v dx7r(x) z 2 D N (x) exists, 
that \zb,n{x)\ ~ 7i (x) everywhere and, lastly, that the integral of tt(x) be known. 

As a first step towards the evaluation of Z D N , we compactify zd^ onto the surface of 
the (i-dimensional unit hypersphere S d . To do this, we introduce polar coordinates, which 
we succinctly denote by (xi, . . . , Xd) = (Qd, R)- The action S 



S 



1, 



D 



is homogeneous in x: 

s(n d ,R) = s(n d ,i)xR 4 

as is the Pfaffian (k = (D — 2)(N 2 — 1)) 

v D , N (n d ,R) = v D , N (n d ,i) x R k 

The R— integration can now be performed exactly, and we find 

fVtt d z D>N (tt d ) 



Z 



D,N 



with 



r((iv 2 - 1) 



2 J S{u d ,l) 



(37) 

(38) 
(39) 

(40) 
(41) 



Without being obligatory, the compactification has nevertheless considerably smoothed the 
integrand, whose i?-dependence was extremely strong. Eq . (|4lD means that we can evaluate 
the integral eq.(^) as the average value of z^^i^d) over randomly chosen points Q d on the 
unit sphere S d . Such points can easily be generated from Gaussian random numbers @ 
Xi, i — 1, . . . , d, projected onto S d . 

We now face an a priori formidable technical difficulty: One needs an efficient algorithm 
for evaluating ZDN^d)- However, it is well-known that Pfaffians belong to the enfants 
terribles of numerical analysis. Direct evaluation of eq.(|5|) is impossible. This has to be 
compared to the case of determinants, where very powerful algorithms are available. There is 
a direct, well-known relationship between Pfaffians and determinant since one has (PLM) 2 = 
det Ai. Taking the square-root, we can always find the Pfaffian up to a sign. Unfortunately, 
we have found that for N > 2 the Pfaffiansf] are not positive semi-definite! As discussed 
in detail earlier, the direct computation of the Pfaffian is possible in generality only for 
D — 4,6. For D = 10 more ingenuity is required, except for SU(2), where a convenient 



representation is available, cf eq. (|T9|) . Using 5*0(10) invariance, we have also found a way 
around the obstacle in the case of D — 10, SU(3): cf eqs.(0),(J3]). D = 10 and iV > 4 
require new ideas, which will not be treated in this paper. 

After computation of z B jv(fi rf ), the evaluation of Zd,n with the simple "direct sampling" 
algorithm is straightforward, as soon as the variance of the integral exists and is sufficiently 
small (cf eq . (|35|) ) . This latter condition is far from innocent, as can be seen from the data 
in table 1, where we present our results for SU(2). We actually give the mean in units of 
V8n, cf eq.(||). 



3 An exception seems to be dimension D 
semi-definite for all N. 



4, where we conjecture the Pfaffian V^n to be positive 
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Table 1: Data for SU(2), direct evaluation of integral 



D 


samples 


Mean 


error 


4 


4 x 10 12 


0.242 


oo 


6 


2 x 10 12 


0.250 


< 0.001 


10 


1 x 10 8 


1.250 


< 0.001 



Indeed, we have been unable to compute the D = 4 integral, because the variance of the 
integral does not exist in that case. Let us mention in passing that the (absolute) divergence 
of any integral / dxf(x) is best checked not by direct sampling as in eqs.(j^),(f40l), but by 
performing a separate Metropolis random walk, to be discussed below, with 7r = |/|. For 
the case of / = z 2 D=iN=2 -, the Markov chain of this random walk quickly gets stuck in one 
valley, which proves the divergence of the second moment for all intents and purposes. For 
the integral itself the random walk finds its way back out of the valley (since the integral 
is finite), but it cannot be calculated by the present method. This extremely powerful 
method to ascertain the absolute convergence of smooth high- dimensional integrals plays 
no role in usual statistical physics applications, where the weight function tt(x) (Boltzmann 
weight or density matrix for classical and quantum applications, respectively) always assures 
convergence. With this simple method, clearcut answers can be obtained both for the 
absolute convergence of the integral itself and for the variance. In cases of doubt, the 
convergence check has to be done first since the relative weight of rare but strong fluctuations 
is difficult to assess by direct sampling. 

Direct sampling of the integral, while successful for D = 6, 10 for £77(2), as shown, has 
proven unfeasible for higher gauge groups. We thus need to consider importance sampling. 
The canonical way to proceed would be to use an analytically tractable approximation 
~ \ z d,n\ with known integral V n . Such functions are unfortunately yet unknown to us. 
Short of this, we will present in the remainder of this paper an efficient algorithm which 
allows us to compute the ratio Z D > ^ N / Z D <^ N for any of the pairs 7J < < D > of interest to us. 
The idea of the method is that for any configuration of matrices X\, . . . ,Xd> describing a 
"valley" configuration with large weight, we will be able to find subensembles oiD < matrices 
(among the above) describing a configuration with an important contribution to Zd<,n- 

To implement the idea, we first note that one condition in the above eq. ( f^Of) was rather too 
restrictive. Because of our particular way of compactifying the integral, we may generate 
random points ild< on the surface of the sphere S d< not only from projected Gaussian 
variables in ]R d< , but also from any other volume isotropic with respect to £^< . Surprisingly, 
one perfectly correct way to generate random f2^< consists in taking d < components of 
random vectors on £^> . Generalizing further, we write 

7r(n Xl ,...,x D> ) = -T,^D<,N^x n ,...,x tD< ) (42) 
v M 

here a = . . . , id<) denotes a combination of D < indices out of D > , and {a} the sum 
over all the combinations. The total number v of possible combinations is 

-(E) 
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L x, 



D< 



denote 



In an effort to be completely explicit, we specify that £lx x ,...,x D> and 
points on the surface of Sa> and S d < , respectively. 

The important point is now that the average value of n(Qxi,...,x D> ) is equal to Z D <^ N . 
We have just transformed Zd<,n into an integral of the same dimension as Zd>,n- The 
last step of our procedure consists in using the same set of X\, . . . , Xd> for a simultaneous 
evaluation of both integrands: 



Z d >,n _ IVQ,f>ir(Qxi,...,x D> ) x z D > ]N {n Xl ,„.,x D> ) /^{^x u ...,x D> ) 



N 



x ly ..,x 



(44) 



In practice, we generate points X 1 , . . . , X D > according to ir(£lxi,...,x D> ) with the Metropolis 
algorithm: starting from an initial configuration on the unit sphere Qd> = ■ ■ ■ , x d >), 
we propose a simple local change of two of the coordinates, which conserves the norm 
of the total vector, to arrive at Q d >- The proposed move is accepted with probability 
min(l, 7i(Q' d> )/ir(Qd>)) (for an elementary introduction cf, eg, ||). We also have to cure the 
sign problem which is completely marginal for D = 4, 6, 10. In fact, since 7r in eq.(f42|) is not 
always positive, we write 



/ VVld 7r x z/tt 





71 


x z/\ir\ 




7T 


sgn(vr) 



(45) 



| vi" | is now a positive weight, which can be used for Monte Carlo sampling M. 

As a result, we obtain the ratio between two different partition functions. We now present 
our results on the ratios of the partition functions, both for SU(2), which serves as a test, 
and for SU(3). Obtaining these data required a few days on ten workstations at one of our 
institutes (AEI Potsdam). 

Table 2: Data for SU(2), ratio of partition functions 



D>/D< 


samples 


Mean 


error 


6/4 
10/ 4 
10/ 6 


6.0 x 10 7 
4.0 x 10 7 
3.8 x 10 7 


1.00 
5.00 
5.01 


<0.01 
0.01 
0.01 



Table 3: Data for SU(3), ratio of partition functions 



D>/D< 


samples 


Mean 


error 


6/ 4 


4.7 x 10 7 


0.99 


0.01 


10/ 4 


6.0 x 10 7 


9.95 


0.09 


10/ 6 


4.0 x 10 7 


9.98 


0.02 



For SU(2), the exact values for the three ratios 6/4, 10/4, 10/6 (cf eq.(p3|)) are, re- 
spectively, 1,5,5, which are reproduced. For S77(3), the numerical computation leads to a 
novel, and very strong result: Here the ratios in question are to a surprisingly good precision 
1, 10, 10, which beautifully agrees with (^) for N = 3. 

Unfortunately, the price to pay for coming up with a powerful algorithm was that we 
had to give up computing the integrals themselves. We have not attempted to test our 
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conjecture eq.fl27|) for the group factor J 7 ^ which appears in the proposed form (|26|) of the 
partition functions Z^,^. A direct computation of the integrals themselves belongs to the 
same class of Monte Carlo problems as free energy calculations, which are notoriously dif- 
ficult. As mentioned before, it seems preferable to find a good analytical approximation 
to the integrands with known integrals. On the other hand, we believe that the difficulties 
in the calculation of the D = 10 Pfaffians can be overcome, and that our ratio test of the 
extended Green-Gutperle conjecture can be extended to gauge groups with N > A. In any 
case, and independent on whether the extended conjecture will turn out to be correct, we 
feel we have already demonstrated the usefulness of numerical calculations in getting reliable 
information on these fascinating and intricate supersymmetric models. 



Note added: We are now able to numerically compute the integrals for SU(3), S77(4) 
and SU(5) directly |I| and to confirm the corrected formula eq . (|26|) . Meanwhile we have 
also received the impressive paper flTT| , where under certain assumptions (i.e. contour pre- 
scriptions) the result for arbitrary N (but without working out J 7 ^) is arrived at by analytical 
methods. 
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